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Introduction 

This preprint represents first, preparatory, stage of numerical investigation of magneto- 
static spin waves' chaos in ferrite films realized between 2001 and 2003 under particular 
support of Multimanetic Solutions Ltd. Its main contents will be submitted in separate 
preprint, as Part III of the manuscript whose Part I was devoted to theory of magnetostatic 
waves and already presented by our arXiv preprint 1204.0200 . 

Our primary purposes here, in Part II, and in the next Part III were: (i) to test numeric 
algorithms based on spatial discretization of film's volume; (ii) to estimate an extent to what 
the theory developed for idealized infinite films is applicable to real finite-size (and may be 
defective) films; (iii) to examine theoretical concepts, - e.g. "quasi-local magnetic energy 
density" (see below), - which have no unambiguous theoretical definition (because of strong 
non-locality od dipole interactions or by other reasons) but can be useful for description 
of non-linear magnetic-wave phenomena; (iv) to visually indicate mechanisms and forms of 
magnetostatic-wave chaos, (v) to see what of them are most appropriate for practical use, 
and numerically investigare possibilities of control and synchronization of this chaos. 

6. NUMERICAL METHODS AND NUMERICAL SIMULATIONS OF EXTER- 
NALLY DRIVEN FILMS 

The modern theory of nonlinear wave processes in ferromagnets and ferrites is not devel- 
oped to an extent sufficient for fruitful applications to so complex phenomenon as magnetic 
chaos. By this reason, numerical simulations must be in anyway useful, since they 

(i) serve as powerful "microscope" to watch for magnetization dynamics at its natural 
temporal scales, from the one tenth part of nanosecond up to tens microseconds, and in this 
way allow to 

(ii) verify existing theoretical models, find prompts for improving analytical theory, and 

(iii) obtain concrete practically acceptable estimates, conclusions and predictions. 
Although many aspects on nonlinear wave dynamics may be modelled in terms of popular 

NSE (nonhnear Shrodinger equation, see Sec. 5), it does not approach for above purposes, 
because it contains no natural amplitude limit for waves and solitons. Indeed, in Eq.ll 
itself and in its solitonic solutions (5.25) and (5.26) the amplitudes, |^| and A , may be 
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arbitrary large, while physically, according to relations (5.9), they can not exceed a level 
of order of unit. This means that at certain critical time moments of chaotic dynamics 
other higher-order non-linearities play significant role. Besides, NSE neglects odd powers of 
non-linearity and related parametric processes. 

Therefore, numerical simulation algorithms should be based on complete Landau-Lifshitz- 
Gilbert equation (2.1). At this approach, trivial normalization of spin length to unit ensures 
taking into account all orders of non-linearity. 

6.1. NUMERICAL ALGORITHM. 

To solve Eq.2.1, the classical third-order Runge-Kutta algorithm was used, in its adaptive 
version with time step being automatically chosen as large as possible at fixed precision. 
Nevertheless, running of Eq.2.1 appears not a fast procedure. The matter is that non-local 
dipole interaction should be calculated by means of Fourier transform, which even in its 
most fast form takes much greater time than, say, calculation of gradients and Laplasians. 

Moreover, if the transform was located within the sample volume then the non-local 
interaction would be represented by a function of two spatial points, ri and r2 instead 
of only their difference ri — r2 , and Fourier transform could not be applied. Therefore, 
when discretizing the sample volume into N^x NyX cells ("spins"), it is necessary 

to consider at least N'^xN'yX N'^ cells, - with N'^^y^^ > 2Nx,y,z — 1 , - in order to correctly 
represent the dipole interaction kernel in the (rl — r2) -space (fortunately, N!^^y^^ — ^N^^y^z — 
1 or N'^^y^^ — 2Nx,y,z are always sufficient numbers). The award for this complication is 
that all the finite-size and boundary effects are taken into account. 

The algorithm was tested by modeling magnetization of small ferromagnetic slabs, with 
sizes of order of tens tq (exchange interaction radius). It was possible to observe (i) forma- 
tion of magnetic vortices and domains at sufficiently weak field, (ii) their death at moderate 
fields, with residual demagnetization at former domain boundaries, (iii) almost uniform mag- 
netization at strong fields with significant demagnetization at sample boundaries only, (iv) 
hysteretic effects, characteristic hysteresis curves on H — B -plane, re-magnetization chaos 
and noise at periodically varying field (chaotic hysteresis). 

6.2. CONCRETIZATION OF NUMERICAL TASKS. 

We will be most interested in comparatively large-scale phenomena in YIG films whose 
typical size is about 7 mm x2 mm xlO micron. Besides, we want to apply such 
kind of film magnetization which ensures maximum MW frequencies at minimum value of 



3 



magnetizing external field Hq . This requirement is satisfied if (i) film is tangentially 
magnetized and (ii) the surface magnetostatic waves (MSW) are explored which propagate 
perpendicularly to static magnetization (see Sees. 3-4). Practically suitable wave length, A 
, of these MSW lies in the interval 0.1 ^ 1 mm , that is X/D > 10 { D is film thickness), 
while suitable cross sizes of MSW inductors (antennas) are comparable with D . 
6.3. ASSUMPTIONS AND SIMPLIFICATIONS. 

If we took the ro (exchange radius ~ 5 • 10"^ cm) be the scale of spatial discretization, 
then (at above mentioned film sizes) A^^x NyX would be ~ 10^^ . Hence, literal 

simulation of real samples is impossible. 

The natural possibility to simplify numerical problems arises from the fact that surface 
MSW with X/D > 10 are almost uniform across film's thickness (i.e. in z -direction). 
Therefore, we can use the averaging over thickness. Then the latter becomes basic spatial 
scale, and we come to two-dimensional discretization lattice, with — 1 and x Ny 
determined by the film length to thickness and wide to thickness ratios. 

In this approach, we inevitably neglect bulk MSW modes with N > which are 
essentially non- uniform in respect to thickness (see Part I, Sec. 4). These non- uniform bulk 
waves are not directly excitable by thick antennae, but they can be generated from surface 
MSW by means of nonlinear fourth-order G-P-processes or third-order P-processes (Sec.5) 
and thus infiuence their dynamics. Nevertheless, our formally crucial simplification has 
physical grounds as follow. 

In fact, there are similar fourth-order interactions between surface MSW themselves. We 
can suppose that just these interactions are dominating in surface MSW dynamics because 
realize in resonant way, while interaction with most of non-uniform bulk modes is far from 
resonances. Indeed, their frequencies lie below uniform precession frequency, Uu ■ while 
frequencies of surface MSW are higher than ouu ■ If comparing the dispersion law (3.38) 
for DE-waves and Eqs.4.28 and 4.29 for dispersion of bulk modes, one can verify that their 
frequencies are close at 

only. In reality, so short-wave modes must be damped certainly stronger than long surface 
waves. At the same time, the latter effectively influence one upon another around (in k- 
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space) equi- frequency curves (shown in Fig.4a, see Sec.4). By these reasons, we expect that 
fourth-order interaction with non- uniform modes is relatively insignificant. The same can 
be argued in respect to third-order parametric by surface waves (all the more, at liJu/2 < 
cji , i.e. at Hq > 4:1t/3 , this process is forbidden at all). 

At present this is (i.e. was ten years ago) unavoidable simplification, since at > 2A'" 
, that is at 80 , again numerical simulation would be unrealistic even with Pentium- 

IV in our order. Rigorous analysis of role of non-uniform (multi-layered) MW in nonlinear 
long-wave dynamics is the interesting task for future. Now, we yet are forced to deal with 
numerical model rather than with literal numerical simulation. 

6.4. NUMERICAL STRATEGY. 

After the simplification, film thickness takes the role of length unit. Moreover, the size of 
discretization lattice, in units of D , can be chosen nxnxl , which allows to simulate 
the film area nN^Dx nNyD . The permissible additional roughening, n , depends on 
characteristic minimum wave length in a situation under analysis. To speed up calculations, 
we may choose a greater n but then lower it if necessary. 

In fact, the values 1 < n < 8 , 20 < iVj^ =< 170 , and 70 < iV^ < 256 were 
used [ten years ago, while now it is possible to take greater , A^^^ and simultaneously 
A^; > 1 , i.e. fractional n]. Of course, the dipole interaction of discrete cells (effective 
spins) was calculated with taking into account their shape as determined by n . A choice 
of definite "magic" and Ny numbers ensured most fast FFT on a N'^ x N'y lattice 
(with N'^^y > 2Nx,y — 1 ). Nevertheless, typically from 1 to 5 real time seconds were elapsed 
per one period of spin precession, ~ 0.3 ns , since several FFT's and inverse FFT's should 
be performed at each time step. 

At given uniform magnetizing field Hq , firstly static magnetization pattern was cal- 
culated and conserved in memory, then serving as initial ground state for wave and soliton 
structures caused by time- varying currents. 

The anisotropy is what can be simply taken into account with no numerical problems. 
But it evolves several parameters at once. At present stage, we want to obtain numerical 
"reference point data" with minimum amount of free parameters and therefore intentionally 
omit anisotropy. 

6.5. NUMERICAL FRICTION. 

It is well known that the time discretization when numerically solving differential dynamic 
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equations inevitably results in more or less effective friction (or may be negative one) . In our 
case this artifact also takes place leading to energy relaxation even if the friction coefficient 
is put on be zero, 7 = 0. Interestingly, this excess numerical relaxation excellently obeys 
exponential law and therefore works as increase of 7 , 7 — )■ 7e// = 7 + 7rmm • 

The value of '^inum depends on mean time step. Typically, the latter was between one 
thirtieth and one twentieth part of the precession period resulting in ^num ~ 0.0004 . This 
friction is just suitable to simulate good but not best samples. However, it could be made 
lower than 0.0001 if decrease mean time step to about one fortieth part of the period. 
Below, the designation 7 will stand for 7e// . 

6.6. MSW EXCITATION BY WIRES AND LOOPS. 

First of all, excitation of weak magnetic waves in small-area film by wire and loop in- 
ductors was numerically watched for. At present, we confined ourselves by inductors with 
round cross-section and radius greater than D , oriented along y -axis in parallel to mag- 
netizing field. Corresponding current induced field, h{r,t) — {hx{x,y,t),0, hz{x,y,t)} 
, was calculated from usual magnetostatics formulas. For relations between physical and 
dimensionless time and frequency units, see Sec. 3. 

Examples of such the simulations are illustrated by Figs.6a-c. Two rather obvious con- 
clusions do follow from these pictures. 

(i) At given microwave frequency of linear wire or loop current, not a single plane wave is 
induced (as it would be in infinite-size film), but a spectrum of waves with different length, 
including ones with non-zero y -component of wave vector (notice that at Hq — 3 the 
uniform precession frequency cuu ~ 6.83 ). In fact, in Fig. 6a and Fig.6b we observe 
eigen-modes (or compositions of nearly degenerated eigen-modes) of small (finite-size) film. 
Comparison of these figures shows that, naturally, loop inductor ensures better wave selection 
and simpler magnetization pattern. 

(ii) At the same time, the spatial Fourier spectrum of excited pattern can be intelligently 
interpreted in terms of infinite film theory (Sec.4), as illustrated by contour plot in Fig. 6b. 
In this plot, a number of lines surrounding some point of k -plane indicates its contribution 
to summary picture. Clearly, spectrum maxima well agree with equi-frequency curves in 
k -plane of infinite film (Fig.4a in Sec.4). in spite of not large film length to wavelength 
ratio ( ~ 5 ). Hence, so visible characteristic rhombic structures in Figs.6a-b directly reflect 
characteristic slopes of equi-frequency curve responding to the excitation frequency. 
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The Fig. 6c justifies that particular eigen-modes and eigenfrequencies of even small-size 
film may be quantitatively close to waves in infinite film. One particular mode is selected 
by equating pump frequency to that of plane surface (Damon-Eshbach) wave with ~ 
±27r/2/ and ky = which would be generated by the same loop in infinite film. We see 
that result is almost plane wave too, i.e. indeed resonance takes place. Nevertheless, this 
mode contains some contribution from plane waves with kx ~ ±3 • 27r/2Z and ky ^ 
which possess the same eigenfrequency and occur resonantly excitable by the same loop. 

6.7. ROLE OF FILM EDGES. 

In the top of Fig. 7, static distributions of the internal field, Wq , and magnetization in 
small-area film are shown, at moderate vahic of external field, Hq = SM^ . Clearly, in most 
part of the film practically uniform magnetization realizes, with Wq ~ Hq . Substantial 
demagnetization takes place at narrow strips only which adjoin film edges perpendicular to 
external field and have width ~ D (as it was stated in Sec.3). 

Since the internal field is lowered at these demagnetized strips, a local spin precession 
frequency there also is lowered. Therefore, usually these strips take almost none part in 
shaping and propagation of waves, as if spins were partially pinned there. These statement 
can be illustrated by Figs.6a-c. 

In principle, specific edge waves can be excited in the demagnetized regions. However, 
this is rather exotic phenomenon, and it was not a case in practically all of our numerical 
simulations. 

6.8. A SOLITON FED UP BY WEAK CONSTANT PARAMETRIC PUMP. 

One more exotic phenomenon is illustrated at bottom of Fig. 7, concretely, very small- 
amplitude soliton (spatially localized wave packet) created and then supported by extremely 
weak parametric pump. The loop current parallel to external field (i.e. to y -axis) induces 
field with amplitude of its x -component ~ 10~^ (i.e ~ 0.01 Oe in real units) and 
frequency uue — 14.5 which is far out off total MSW frequency band (with upper bound 
— Hq + 27: 9.3 , see Sec.4). As the result, long weak envelope soliton is formed 

whose carrying frequency equals to half of tUe , with magnitude ||S'_La:|| ~ 10~^ and 
width ~ 150D (at film length = 432D ). Interestingly, the two latter quantities 
approximately satisfy the relation between amplitude and width of bright solitons which 
follows from Eqs.5.15 and 5.25. Although bright solitons moving perpendicular to external 
field are formally forbidden in infinite film, this simulation shows that similar objects are 
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permitted for finite-size film. Besides, we detect that parametric excitation allows to create 
small- amplitude soliton avoiding the restriction (5.23). 

Further behavior of this object is even more intriguing. It oscillates between film ends 
(sec Fig. 7) undergoing some decay after reflection from the distant end but more or less am- 
plification while reflecting from the edge where inductor takes place. The dot line separates 
period of strong amplification of the soliton due to occasionally "good" relation between 
phases of its carrier and pump. Dimensionless time and frequency in Fig. 7 are real ones as 
expressed in units of tq and /o , respectively (see Sec. 2 and 3). 

From the point of view of Eq.5.1 the parametric process under discussion must be de- 
scribed by quadratic term {Sq, h) [Sq, S±] in second row. Interestingly, in view of hy = 
we should conclude that pump, (^o, h) , mostly acts at the demagnetized edges. Among 
our numerical collection, this is exclusive example when film edges play a key role. 

Below, we will deal with non-parametric pump whose frequency belongs to the MSW 
band. Like here, in all forthcoming examples inductors are parallel to external magnetizing 
field. 

6.9. CREATION OF SOLITONS BY NON-PARAMETRIC PULSE PUMP. 

The Fig. 8a shows the consequence of intensive radio-impulse of current passed through 
wire inductor at one of ends of relatively large-area film. The impulse duration was about 
twenty periods of the carry frequency, oUe — 7.5 . At a distance from the inductor, the in- 
duced magnetization precession impulse is deformed. If its initial amplitude was sufficiently 
large then further it breaks into a chain of pulses with almost zero dips between highest of 
them. The picture in Fig. 8a can be quahfied as formation of gray solitons inside finite-length 
wave packet. 

The critical breaking level of amplitude is just its maximum after breaking, in this example 
\\S±x\\ ~ 0.15 . This observation is in reasonable agreement with the estimate of this level 
which follows from Eqs.5.9, 5.10 and 5.23, 

\\S±^\\ ^ A„ii„v^ = Vrp/|x| « 0.1 (2) 

(at friction 7 = 0.0005 what took place) . In all the below discussed numerical simulations, 
close threshold values, 0.1 < ||'S'_La;|| < 0.2 , mark transition to brightly expressed non-hnear 
effects and to magnetic chaos. 

6.10. CHAOS UNDER UNIFORM RESONANT MICROWAVE FIELD. 
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In most of previously reported experiments on magnetic chaos, the latter was excited by 
nearly uniform microwave magnetic field, h{r, t) , either using parallel parametrical pump 
when h\\So and excitation frequency ~ 2a;„ or through perpendicular ferromagnetic 
resonance (FMR) when /iXS'o and c^e ~ • Consider the second variant since it is 
more close to chaotic auto-generation of MSW to be under our interest. 

The Fig. 8b demonstrates modeling of nonlinear FMR in moderate-area film at Hq — 3 
and 7 = 0.0007 , under uniform microwave field parallel to x -axis. Due to finite-size 
effects (edge demagnetization), factual (numerically found) uniform precession frequency, 
u!u ~ 6.785 (at Hq — 3 ), is shghtly lower then theoretical value for infinite film, 
^uoo ~ 6.83 . Taking u^, sufficiently close to Uu , it is possible to obtain strong response 
to weak perpendicular field h ~ COOIM^ ~ 0.1 Oe . Naturally, the response is indicated 
by amplitude of uniform component of S± , i.e. {S±) = JyS±dr/V where V stands 
for film volume. 

After sharp two times increase of pump we observe increase of {S±) ' amplitude which 
monotonically tends to nearly two times larger value. However, next increases of h by 
the same step result in smaller and non-monotone response. At h = 0.007 , the response 
transforms into periodic oscillations. At last, when h = 0.008 , these oscillations turn into 
chaotic one, at time moment marked as "burst" on Fig. 8b. More careful repetition of this 
process (with smaller step) allows to notice at least one or two period-doubling bifurcations 
of regular oscillations (as well know in theory of dissipative chaos [1]). 

Interestingly, to reach the chaos, well satisfied resonance condition oue ~ ouu is quite 
necessary. For instance, at LOe = ^uoo (i-e. at frequency deviation less than 1% ) 
the only result of even very intensive pump, h ^ 1 , is strongly nonlinear but regularly 
oscillating long-wave structure. 

Other characteristic observation is the hysteresis of chaos: if h is lowered from 0.008 
than chaotic regime remains at least down to h — 0.005 . 

6.11. SHORT-WAVE EXPLOSION AND TRANSITION TO CHAOS. 

Usually, experimental magnetic chaos is analyzed in terms of particular wave modes, that 
is in momentum space (see, for instance, [2,3]). In our simulation, we can view also how it 
looks in real space, and watch for spatial-temporal picture of transition from regular motion 
to chaos. Characteristic scenario of this transition is illustrated in Fig.9. 

When FMR is still regular, rather smooth magnetization pattern takes place with one or 
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three maximums of oscillations and precise mirror symmetry (with respect to middle lines 
of rectangular film area). However, the closer is the transition the higher and narrower is 
the central maximum. This means that spectrum of excited MSW becomes more and more 
wide, but still coherent, in the sense that all the waves are mutually connected by some 
rigid phase relations. At critical "burst" time moment the central maximium collapses into 
peak very narrow in x -direction and rather flat (elonged) in y -direction (i.e. along static 
magnetization). Then this peak blows up giving freedom to to the short waves. The latter 
incoherently scatter in all the directions in their turn giving rise to complicated (chaotic) 
magnetization pattern. Beginning of this unstable explosive stage is shown at bottom of 
Fiq.9. Notice that corresponding magnitude of S^-j^ confirms estimate (2). 

Important sign of transition to chaos is violation of the mirror symmetry. The symmetry 
with respect to 180° rotation only remains after the explosion. Clearly, this symmetry is 
invoked by that of the static -S'o pattern (see Fig. 7). 

At later time, most short of the explosively induced waves decay. In further stationary 
chaotic regime, magnetization picture more or less restores both smoothness and mirror 
symmetry. But, naturally, increase of pump rises both short-wave contents and asymmetry. 
Moreover, under sufficiently intensive pump new similar explosions (bursts) are repeated 
from time to time, serving as "discharges" of excessive energy accumulated by long waves. 
This may be called strong chaos. 

6.12. EXCESS ENERGY AND POWER ABSORPTION. 

The top of Fig. 9 demonstrates typical chaotic behavior of excess film energy, E , and 
of power absorption by film, P , both related to unit volume (discretization cell). Here 
and below, "excess energy" (or simply "energy" ) will term increase of magnetic energy due 
to the excitation, ^j, , but excluding direct S±_ 's interaction with pump field (i.e. 
except —J {h,S±) dr /V , see Sec. 2). The power absorption, P , describes energy flow 
into film which is spent for both E and dissipation in the film interior. Hence, in general 



with Prfis being dissipated power per unit volume. Symbol designates time 

averaging with respect to spin precession. These equations directly follow from the basic 
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Eq.2.1. 

It is natural to expect that 
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Pdis^2r{E-Eo) , (4) 

at some Eq , and F being previously considered dissipation rate of magnetization. At 
7 = 0.0007 first of Eqs.5.10 gives 2r 0.013 , then energy relaxation time is ~ l/2r fa 
80 . Indeed, approximately such the time scale can be viewed at plot (A) in Fig.9. Further, 
let us put on {S±x) — \\Sx\\sin[uJet — 0] where \\Sx\\ and designate amplitude and 
phase, respectively, of the uniform component of spin precession. Then 

P = h (sm{ujet) ^ ^^^'^^ ^ p=i^uJeh\\S^\\sm 4> (5) 

At ll^'j:!! 0.1 (as prompted by Fig.Sb), /i = 0.008 and ?s 7r/2 (which means good 
resonance) this relation results in P ~ 2.5 • 10~^ , in agreement with plot (B) in Fig.9. 
Hence, rough average characteristics of chaotic variables are easy explainable. 

Much harder task is to explain chaotic deviations from average values. For instance, in 
plots (A) and (B) energy and power nearly follow one another. It is clear: the energy comes 
from power absorption, but the latter depends on the phase which in its turn must be 
sensitive to energy, because of non-isochronity of spin precession (see Sec.5). But without 
an adequate dynamical model for this connection we can not estimate details of chaotic time 
series. 

6.13. INSTANT FREQUENCY. 

Directly, our numerical algorithm produces fast oscillating time series. To extract from 
them relatively slow time- varying ("instant") amplitudes and phases (like \\Sx\\ and 
above), and corresponding instant frequencies (e.g. uie — d(f)/dt ), there are two ways. 
One is to build analytical signals by means of (discrete) Gilbert transformation. Another 

way is to indicate maximums, minimums and zero-crossings of oscillating variable. If t„ , 
n = ...,—1,0,1,... , are estimates of time moments when zero-crossing docs occur, then 
instant frequency at t ^ tn can be determined as uiin — 2k'n/{tn+k — tn-k) , where 
A; > 1 . 

If disposing this quantity,the phase can be restored by discrete numerical integration. 
This method is more comfortable and fast than Hilbert transformation, but ensures not 
worse (usually better) accuracy what was confirmed by special tests. 

6.14. QUASI-LOCAL ENERGY DENSITY. 
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Because of relations (2.9) and (4.1), the excess energy (per unit volume) can be repre- 
sented in the form 



E, 



nonloc : 



E, 



loc 



dr 

V 



= Ho{l - Sy) + 2ttS[ 



(6) 



Here Enonioc is contribution from non-local (non-singular) part of dipole interaction in 
plate geometry (Sec.4), while Eioc consists of its local (singular) part and also local first 
term of Eq.2.9. 

For uniform precession (even let large- amplitude and strongly non-linear), Enonioc 
vanishes, hence, the sum cioc — {h{t), S) plays the role of Hamiltonian of the average spin. 
In particular, in autonomous regime of uniform precession, at h — , eioc becomes 
integral of motion. 

Importantly, Enonioc ma-y be negligible as compared with Eioc in non-uniform chaotic 
case too. Indeed, according to equations (4.1-2), (4.6) and (6), the non-local contribution 

can be estimated as 
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dk 



with Sx,z{k) denoting spatial Fourier transform of magnetization. If long MSW (with 
D\k\ ) are dominating in magnetization pattern, then the excess energy is well char- 

acterizable by Eioc ■ As the consequence, the quantity eioc behaves like local integral of 
motion and thus can be termed quasi-local energy density. 
6.15. WEAKLY CHAOTIC FMR. 

Seemingly weak chaos reahzes aX h = 0.005 . The Fig. 10 shows evidence for (i) good 
correlation between the energy and dissipation (with Eq = Q and F well related to 
actual friction coefficient 7 ), and, at plot (D), (ii) rigid correlation between energy and 
average longitudinal component of magnetization, (^S\\) . Here (..) denotes space-time 
averaging. Notice that in most part of film area S\\ — {So, -S") fa Sy . 

The latter correlation just gives the evidence that non-local energy contribution is rela- 
tively small. Indeed, because of relations 
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11^.11 Vp^ , <^.^> + <^.^>-2(l-(5|l», 
where p is characteristic eccentricity, Eioc can be expressed as 

Eioc-^{h, + ^^{1-{S\^) (8) 

Hence, at Enonioc <^ -^Zoc total excess energy also well reduces to (^n) . If we equated 
p to small-amplitude uniform precession eccentricity given by Eq.5.10, then at Hq = 3 
the Eq.8 would yield dE/d{S\\) ?s dEi^Jd{S{) fa -5 . This value is lower by 10% 
than the slope ( ^ —5.6 ) at Fig. 10. D. The difference can be explained if take into account 
that increase in spatial non-uniformity of precession result in decrease of its eccentricity. 

Plot (G) in Fig. 10 presents spectrum of power absorption (absolute value of Fourier 
transform of P{t) ). At A/, ^ 140 Oe (as for YIG), the dimensionless frequency unit 
corresponds to ~ 390 MHz, hence, the dominating frequency in P[t) 's spectrum is 

20 MHz. Chaos in instant frequency of precession, as shown in plot (B), and thus in 
its phase is characterized by significantly wider frequency band (in part contained by plot 
(C) in Fig. 9). Nevertheless, plot (F) in FiglO visualizes anti-correlation between instant 
frequency and energy which corresponds to negative sign of non-isochronity (see Sec. 5). 

For fractal dimension (Sec. 5. 9) of these data it was found that dcor < 3 (see below). 
This implies that chaos is governed by three relevant variables only, and therefore chaotic 
attractor could be represented in 3-dimensional space. It seems doubtless that E , P 
and \\Sx\\ are relevant variables (due to Eq.5, P and determine also the phase 

(f) ). Plot (C) in Fig.lO gives show of the attractor in these coordinates. 

6.16. ANALYSIS OF CORRELATION DIMENSION. 

Let us discuss practical calculation of the correlation dimension, dcor (Sec. 5. 9). At given 
correlation sum, one may estimate dcor by two ways: 

_ c^ln a{R) _ \n[a{R)/n] 

The first of them may be called differential dimension while the second integral. Under 
formal limit N ^ oo these quantities are expected to coincide one with another. But 
real calculation needs in more than N'^ operations, therefore N can not be as large as 
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wanted. Under realistic N , the differential estimate satisfactorily works at moderate 
values of R only, i?min <S R <^ -Rmax (where i?min is minimum of Rij ) . In 
opposite, the integral estimate better works at lower end of this interval. 

Naturally, the integral estimate is less sensitive to finiteness of N , but instead it requires 
to know the coefficient fl . It is rather obvious that the best general assumption about fl 
is that it equals to volume of dear -dimensional unit-radius sphere, i.e. il — 7r'^™''/^/r(l-|- 
dcor/'^) , where F means gamma-function. By special tests we verified that this recipe 
indeed constantly improves precision of dcor estimates. 

Notice that wide class of tests is presented by the Kaplan- Yorke chaotic system [1] . It is 
described by the set of difference (discrete time) equations: 

n 

X,{t+l) = F{X,it)) , j^l...n, r(t + l) = ar(t) + J]/,(X,(t)) , |a|<l, (10) 

where F{X) is some chaotic one-dimensional map (for instance, tent map, F{X) = 1 — 
\2X — 1| ), and fj{X) arc any smooth functions. Dcpcndcntly on n and a , fractal 
dimension of chaotic sequence Y{t) can be equated to arbitrary number (for example, 
for the tent map n — 2 and a = 1/16 lead to dcor — 2.5 while n — 3 and a — 1/64 
to dcor — 3.5 ). 

6.17. FRACTAL DIMENSION OF CHAOTIC FMR. 

Typical example of evaluation of dcor is shown in Fiq.lla (left plot). Here "log of Cell 
Size" in horizontal axis means In (Rmax/R) while vertical axis presents both the estimates 
(9). The dependencies of differential and integral dimensions on In (i?max/-R) are drawn 
by thin and fat lines, respectively. 

We see that plateau in the first of them well coincides with the upper value of the 
second. From above mentioned tests, it is known that just this value should be taken as 
best estimate of dcor , a-nd that the coincidence signify good reliability of this estimate. 
In more complicated case when the plateau differs from upper (most right-hand) value of 
integral dimension, the latter must be preferred. But such the situation testifies that either 
the data (finite chaotic series) are not enough representative or they possess essential multi- 
fractality. 

In case under consideration, fractal dimension of the power absorption time series under 
weakly chaotic FMR can be estimated as dcor ~ 2.4 . 
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To form better representative data, the time separation, r , in d -dimensional embed- 
ding points {x{to + nT), x{tQ + nT + T), ...,x{to + nT + dT)} should be a few times shorter 
than characteristic correlation time of x(t) , while number of points must be sufficient for 
minimum "filling" of all of d dimensions, at least N >2'^ . 

6.18. CHAOTIC FMR IN DEFECTIVE FILM. 

Real films always have more or less amount of defects. The right-hand plot in Fig. 11a 
demonstrates static magnetization by tangential field in film with periodic lattice of defects 
(punctures) which touch about 10% of film area. Naturally, mean internal field, Wq , and 
thus characteristic precession frequencies are lowered by defects. But we found no qualitative 
difference between chaotic FMR in defective film and "good" film. 

The top and left bottom plots in Fig. lib illustrate how the uniform component of spin 
precession behaves at small time scale (of order of Ins ) and at moderate time scale (of 
order of 0.1 //s ). On right hand, contour plot of spatial Fourier transform of S_i_x shows 
that two groups of MSW modes form magnetization pattern, both belonging to the same 
frequency range around uiu , with the short-wave group induced by defect lattice. The levels 
for this plot were chosen specially to highlight short-wave modes. In fact, their contribution 
to energy is of order of a few percents. 

6.19. SYNCHRONIZATION OF CHAOTIC FMR. 

In the work [4] synchronization of magnetic chaos under non-linear FMR in normally 
magnetized YIG film was experimentally realized. First, the signal, Pm(^) , related to 
power absorption was recorded into a memory. Characteristic frequencies of this signal were 
between 0.5 MHz and 10 MHz. The similar actual "slave" signal, Psit) , was compared 
with the recorded "master" signal, and the difference was directed, with some proportionality 
constant, K , to perturb the external magnetizing field, Hq ^ Hq + K{ Pm — Ps) ■ 
At suitable choice of K , after a transient time ~ 10 //s , excellent coincidence between 
Ps{t) and PM{t) was observed. 

The Fig. 12 illustrates the attempt to numerically reproduce such the experiment but 
with tangentially magnetized film. The above mentioned defective film model is explored, 
at Hq = 3 and h = 0.005 . The power absorption per unit volume, P{t) , is taken 
to serve as the control signal. The master signal to be addressed to the feedback, Puif) , 
either equals to P{t) or formed from it by slight time-smoothing (over 5 -^ 30 periods of 
precession) . Typical magnitude of chaotic P{t) 's variations is ~ 5 • 10~^ . The feedback 
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coefficient, K , is changed in the range between —30 and —600 . 

Unfortunately, 10 /is is rather large time for our numerical simulations. Total duration 
of numerical runs was just about 10 /is . Best signs of synchronization were observed at 
K ~ 200 . Hence, magnitude of the bias field modulation was ~ 0.01 , i.e. ~ 1.5 Oe in 
real units. For comparison, in [4] essentially smaller values ~ 0.1 Oe were in action. 

According to Fig. 12, with no doubts synchronization takes place, but its quality is far 
from so excellent as reported in [4]. Possibly, this is due to wider frequency range of the 
control signal in our system, up to ~ 50 MHz (see plot of P{t) 's spectra in Fig. 12) and 
to not long enough duration of the numeric experiment. 
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Conclusion 

The above expounded material leads to conclusions as follow: 

(i) Concepts, formulas and observations of linear theory of magnetostatic spin waves 
(MSW) in infinite films (Sections 2-4 of Part 1, see arXiv preprint 1204.0200), as well as 
that of quasi-linear MSW theory (Sec. 5 in the Part 1), appear quite adequately useful for 
interpretation of numeric simulations of even non-linear and chaotic MSW even in small-size 
films; 

(ii) Most effective way to MSW chaos is via parametric resonance and parametric non- 
linear transformations of MSW; 

(iii) MSW chaos typically is determined by a few relevant variables only, - i.e. character- 
ized by rather low fractal dimension , - and therefore seems allowing its control and more 
or less satisfactory synchronization. 

The two last features will be further investigated in next Part III of this manuscript. 
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concentrating on not externally driven but auto-generated MSW chaos. 
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Fig. 6a. Instant S^^ ^-pattern in small-area in-plane magnetized film under weak 



non-resonant excitation by wire positioned at x/D=9. Both external field and wire are 
parallel to y-axis. Parameters: 1^=3, y=0.003, exciting frequency 0)^=7 . 
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Fig. 6b. (Top) The same as in Fig. 6a but inductor is loop taking lines x/D=9 and x/D=17. 
Both external field and loop are parallel to y-axis. Parameters: F|^=3, y=0.003, (0^=7.3. 

(Bottom) Contour plot of spatial spectrum -pattern. The curve marked by arrows is 

theoretical equi-frequency line co=cOg=7.3 for surface waves (see Fig.4a). 



20 
15 



10 



10 



20 



30 



40 50 

x/D 



60 



70 



80 



20 



15 

N 

^10 



.\\\, 

_\ s . . V 





\ 



^^^^^ 



V I 





_! ^ UA. I 1 



vv^^^ I 

J iJ^ .V I ^ ^ 




10 



20 



30 



40 

x/D, S 



50 



60 



70 



80 



_Lx 




Fig.6c. (Top) Instant S^^^ pattern under nearly resonant excitation by loop current with 
frequency co^=8.26 =cOp^(27tD/2l), l=8D. Other parameters as in Fig. 6b. 

(Middle) Corresponding instant distribution of the vector S^^. 

(Bottom) Contour plot of spatial spectrum of S^^^. The curve marked by arrows is 
equi-frequency line corresponding to co=(o (see Fig.4a). 
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Fig. 7. (Top & Middle) Internal field, W^, and x-component of static magnetization, S^^, in small-area 
film under external field Hq=3 parallel to y-axis. (Bottom) Parametric subharmonics generation under 
excitation by weak loop induced field h^=±O.OOC0(^in(co^t) applied to left edge, at Hq=3, y=0.0005, 
(0 =14.5. Left: magnetization disturbation in the middle of film. Right: its spectrum occurs nearco 12. 



Q 
O 
o 

CM 



0.2 



0.1 



-0.1 



-0.2 



30 





40 



50 



60 



70 



80 



90 



t (time, in units) 

Fig.Sa. Beginning of soliton formation under strong wire-induced field with h^~0.2sin((»^t) 
applied to left edge, at Hq=3, y=0.0005, (0^=7.5. The envelope of Sj^^(x=200D,t) is shown. 
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Fig. 8b. Nearly resonant excitation of in-plane maanetized film by uniform radio-frequency field, 
h^=h sin((Ogt) , at Hq=3, y=0.0007, (0^=6.785. (Top) Transition from stable FMR regime to chaos 

as the field, h, grows from 0.001 to 0.008. (Bottom) Chaotic oscillations of uniform 
magnetization component <Sj^>=j^Sj^(r,t)dr/V2 jn both plots its amplitude is shown. 
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Fig.9. Chaotic ferromagnetic resonance of in-plane magnetized film excited by uniform 
radio-frequency field, h^=0.008 sin(o)gt), at H^^S, y=0.0007, ©^=6.785. (A) Excess energy 

per unit volume. (B) Power absorbtion per unit volume. (C) Spectrum of particular chaotic 
realization of <Sjj^>. (D) Instant S^^ patte^^ before transition to chaos. (E) Characteristic S^^ 

pattern at the moment of transition to chaos (marked "burst" in Fig.Bb). The explosive birth 
of short maanetic waves and mirror svmmetrv violation are clearlv seen. 





Fig. 10. CInaotic FMR in uniform field, ln^=0.005 sin(a)gt), at Hq=3, y=0.0007, (0^=6.785. 
(A) Excess energy. (B) Instant frequency of <Sj^x^*)^ determined by its zero crossings. 
(C) Fragment of chaotic attractor. (D) Correlation between energy, E, and <^ |>. 

(E) Correlation between energy, E, and instant frequency. (F) Correlation between 
dissipated power, P^.^, and energy. (G) Sggctrum of power absorption. 
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Fig.11a. Left: Fractal (correlation) dimension of absorbed power signal, P(t), obtained 
by its embedding into 10-D space. The result is d^^ =2.4 . Right: S^^-component of 

static magnetization in periodically punctured film. 




Fig. 1 1b. (Top) Small piece of uniform component of magnetization excitation, <S^^> (dotted 
lines mark periods of exciting field). (Bottom left) Fragment of <Sj^^>'s envelope in chaotic FMR 
regime of in-plane magnetized film. (Bottogi^right) Contour plot of typical S^^^'s pattern. 
Dotted lines are theoretical equi-frequency lines for co=(o^, i.e. |ky|/|k^|=(47t/HQ)^^^, at Hp=3. 




Fig. 12. Synchronization of cinaotic power absorption signal, Pgi^^gW. at one film undergoing FMR 
by power signal, P^igstgi-W' f''orn another film also subject to chaotic FMR. The films are identical, 
both periodically punctured, tangentially magnetized by field HQ=3Mg and excited by the same 
uniform field, h^=0.005MgSin(cOgt), whose frequency nearly equals to the uniform precession 
frequency, co =5.65. For concreteness, in plots the value M =140 Oe is substituted. 

25 



